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Abstract - We study the effects of random pinning on the Fredrickson-Andersen model on the 
Bethe lattice. We find that the nonergodic transition temperature rises as the fraction of the 
pinned spins increases and the transition line terminates at a critical point. The freezing behavior 
of the spins is analogous to that of a randomly pinned p-spin mean-field spin glass model which has 
been recently reported. The diverging behavior of correlation lengths in the vicinity of the terminal 
critical point is found to be identical to the prediction of the inhomogeneous mode-coupling theory 
at the As singularity point for the glass transition. 


Introduction. — The nature of the glass transition 
is still elusive It is a major challenge to under¬ 

stand universal behavior of the transition, such as the 
non-Arrhenius growth of the viscosity, non-exponential re¬ 
laxation of the correlation functions, and spatially hetero¬ 
geneous dynamics of the fluctuations. The ultimate goal 
is to establish whether the transition is associated with 
any thermodynamic singularity, or it is caused by purely 
kinetic mechanism m- Obvious reasons that hamper 
the progress of our understanding of the glass transition 
are the inhibitedly long time scales required to reach the 
transition point experimentally and the lack of an ideally 
simple finite-dimensional model which exhibits a true glass 
transition, if any. 

Recently, a novel idea to bypass the difficulty to access 
the glass transition temperature by randomly freezing, or 
pinning, a fraction of degrees of freedom of the equili¬ 
brated system has been proposed HHZ]. Cammarota et al 
have analyzed the effects of random pinning on the glass 
transition of the p-spin mean-field spin glass model EC]- 
It was found that both the ideal glass transition tempera¬ 
ture, Tx, and the dynamic transition temperature, Td, rise 
as the fraction of the pinned spins, c, increases. Further¬ 
more, the two transition lines, Tk{c) and Td{c), are found 
to merge and terminate at a finite c. This end point is 
argued to be a critical point whose universality class is 
that of the mean-field random-field Ising model [ 8 ] . In the 
terminology of the mode-coupling theory (MCT), which is 
a mean-field dynamical theory of the glass transition, this 


end point is characterized as the A 3 singularity, where the 
anomalous dynamical scalings, such as the logarithmic re¬ 
laxation dynamics and distinct diverging length scales, are 
predicted [9]. Verification of the ideal glass transition of 
randomly pinned systems by experiments and simulations 
for realistic systems may be crucial to prove (or disprove) 
the very existence of the bona-fide glass transition point at 
finite dimensions nan]- At the same time, it is imminent 
to establish the relationship of thermodynamic scenarios 
such as the random first order transition theory (RFOT) 
with other kinetic scenarios which do not necessarily re¬ 
quire thermodynamic singularities behind the glassy slow 
dynamics pl fT^HTl] . 

In this letter, we argue that the singular behavior of 
the glass transition by random pinning analogous to that 
of the p-spin mean-field spin glass model can be also ob¬ 
served for a purely kinetic model on the Bethe lattices (or 
random graphs). Kinetically constrained models (KCMs) 
are toy lattice or spin models in which Hamiltonian is free 
from the interaction but a non-trivial constraint is im¬ 
posed on the dynamic rule [2j[Tll[T5]. Many KCMs are 
known to display slow dynamics in the collective move¬ 
ment of the particles or spins at low temperatures but they 
are of purely kinetic origin since their thermodynamics are 
trivially ideal-gas like. Their glassy behavior is genetically 
distinct from those of thermodynamic scenarios such as 
RFOT where thermodynamic singularities encoded in the 
free energy landscape play a pivotal role and escort the 
slow dynamics. Indeed, most KCMs show no glass transi- 
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tion at finite temperatures. Notable exceptions are KCMs 
on the Bethe lattices (or random graphs) [T4l[T6]. One of 
the KCMs called the Fredrickson-Andersen model (FAM) 
is known to undergo a nonergodic transition called the 
Bootstrap Percolation (BP) transition at a finite temper¬ 
ature Td on the Bethe lattice UMIIl. It has been demon¬ 
strated that the time evolution of the order parameter of 
the FAM shows dynamical anomaly similar to those pre¬ 
dicted by the MCT for structural glass formers, such as 
the two step relaxation and algebraic divergence of the re¬ 
laxation time near Td [16]. Here, we demonstrate that the 
randomly pinned FAM on the Bethe lattice also shows the 
analogous critical behavior as that of the randomly pinned 
p-spin mean-field spin glass model such as the A 3 higher 
order singularity ism. 

Model. — We consider a FAM on the Bethe lattice. 
In the absence of pinned spins, this model is defined as 
follows. We consider N Ising spins G { —1,+1}, i = 
1, • • • ^N. The Hamiltonian of the system is free from the 
interaction and can be written as 

1 ^ 

( 1 ) 

i=l 

The equilibrium distribution of the up spins is given by 

P 1 i 7 -I /rT-i\ 5 (^) 

1 + exp(-l/T) 

where p varies from 1 /2 to 1, corresponding to T = 00 and 
T = 0, respectively. 

The spin variables represent coarse grained mobility 
fields. ai = +1 represents the immobile region and 
(Ji = — 1 is the mobile region or defect m- The idea 
of the FAM is to introduce a kinetic constraint on the 
time evolution of the spins in such a way that flipping of 
a spin is more difficult when it is surrounded by many up 
spins (immobile regions). To be more specific, the Tth 
spin at each Monte-Carlo step can flip with a transition 
probability 

w{ai -ai) = min{l,e“°'‘/^}, (3) 

only if the number of the nearest spins in the state —1 is 
larger than or equal to / pro] - It is known that the FAM 
on a lattice in finite dimensions does not exhibit the tran¬ 
sition at a finite temperature but it does so on the Bethe 
lattices and random graphs HU- The FAM undergoes the 
nonergodic transition at a finite temperature if the con¬ 
nectivity k, or the number of the nearest neighbours + 1 , 
satisfies k > f > 1 m- 

Now we consider to pin or freeze a fraction c of the spins 
randomly chosen from the N spins which are initially equi¬ 
librated at a temperature T before pinning. Note that 
pinning spins in the equilibrated configuration at T is es¬ 
sential in the following argument. The dynamic rule after 
pinning is basically the same as the bulk system; the i- 
th spin can flip with the probability of eq.(|3|), only if the 



Fig. 1: The Bethe lattice with k — 

number of the nearest spins in the state —1 is larger than 
or equal to / and the i-th spin is not pinned. In the fol¬ 
lowing, we only consider the case of /c = 3 and / = 2 (see 
Fig. Hj), but conclusions for different sets of (/c,/) do not 
change qualitatively as long as /c > / > 1. 

Phase diagram. — Long-time behavior of the ran¬ 
domly pinned FAM on the Bethe lattice starting from the 
equilibrium distribution at the initial time can be eval¬ 
uated analytically invoking the recursive relation on the 
lattice. We shall evaluate the persistent function 0, which 
is the total fraction of the dynamically frozen spins at the 
long-time limit, defined by 

0 = c -h P+ + P-. (4) 

Here the first term c represents the fraction of the pinned 
spins and P± represent the fraction, or the probability, of 
the spins ultimately arrested in the +(—) state among the 
{1 — c)N unpinned spins, respectively. P+ is written as 

P+ = (1 - c)p^ (^ + (5) 

n=0 ^ ' 

where p is given by eq. (|2|) and B is the conditional prob¬ 
ability that a spin is in the state —1 or flipped down to 
the state —1 given that one of the nearest neighbours was 
arrested in the state -hi. B obeys the following expression; 

f /h\ 

P = 1 - p -h (1 - c)p ^ I 1 P^-^(l - BY. (6) 

n=0 

The first term, 1 — p, on the right hand side of eq.([6|) 
represents the equilibrium probability of —1 spins. The 
fact that it is independent of c is the reflection that the 
spins have been randomly pinned from the equilibrium 
distribution. The second term represents the probabil¬ 
ity of the spins originally at the state -hi but eventually 
flipped down to the state —1. Eq. (|6|) is closely related 
to that of the multi-component extension of the FAM and 
associated models in Refs [111120]. Indeed, eq.® can be 
obtained by freezing one of the degree of freedom in the 
equation for the binary system considered in Ref [20] (see 
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Fig. 2: p as a function of B for /c = 3 and / = 2. Three lines 
correspond to c = 0.1, 0.2 and 0.3 from the top to the bottom. 
Inset: Same results for /c = 15 and / = 7, where c = 0.5, 0.6 
and 0.7 from the top to the bottom. 

eq.(3) therein). Likewise, P_ can be written as 



Fig. 3: p dependence of (0(c,p) — c)/(l — c) for several c’s: 
The transition is discontinuous at c < 1/5 and continuous at 
c = 1/5. For c > 1/5, only the crossover is observed. Inset: 
The t-dependence of 0(t) for c = 0.1 and several p’s obtained 
by MC simulation (filled symbols). The solid lines are the 
solutions of the stationary equation, eq. 


f—\ 

P_ = (1 - c)(l - p) ^ f ^ + L P'y 1 - (7) 

n=0 V ^ / 

where B' is the conditional probability that a spin is in 
the state —1 or flipped down to the state —1 given that 
one of the nearest neighbours was arrested in the state — 1 , 
which is given by the solution of 

B' = l-p+{l-c)p Y, ( ( 8 ) 

n=0 


From eqs.Q and 0 , the total fraction of the arrested 
spins, eq.dH, is written as 


</, = C + (1 - c) (B) + (1 - p)^i^, (B') 


(9) 


where we have defined an auxiliary function 


/-I 




n=0 


Q^(l - Bf-^B^. 


( 10 ) 


One observes that B' and 0 can be computed from the 
self-consistent equation for B given by eq. In the limit 
of c = 0, these expressions reduce to those studied in Ref. 
[E]. Eq. ([9]) obtained for the stationary state coincides 
with the long-time limit of the time-dependent persistent 
function 0 (t) evaluated for the equilibrium initial configu¬ 
ration p^fTOfTSlfT^ . We checked this by the Monte-Carlo 
(MC) simulation for several state points (see the inset of 
Fig. El). For = 3 and / = 2 , eq.® can be explicitly 
written as 


B = l-p+{l-c}p{B^ + 3B\l-B)}. (11) 

If c = 0, there is always a trivial solution B = 1 , z.e., 
all spins can flip eventually and the system is ergodic. 
However, when the temperature is lowered from above, or 


p is raised from below from p = 1/2, B discontinuously 
jumps from 1 to = 1/4 at a transition point, pd = 
8/9, and concomitantly (p jumps from 0 to a finite value 
of (pd = 0.673. This is the A 2 transition in the MCT 
terminology [3]. 

If c 0, P 7 ^ 1 always holds and eq. m can be solved 
immediately to obtain 

2 _ ^ 

^ l-(l-c) [P3 + 3P2(1-P)]' 

In Fig. [21 we show p{B, c) as a function of P for several c’s. 
If c 7 ^ 0 but small, P-dependence of p shows a nonmono¬ 
tonic behavior (see c = 0.1 of Fig. |2j). When p is small 
(high temperature) where most of the spins are down (or 
mobile), there is only one solution for P close to 1, where 
a small fraction of the spins freezes in the vicinity of the 
pinned spins and most of the spins can still flip. As we 
increase p, the nontrivial branches at a small P (or a large 
p) appear discontinuously at p = which is the signal 
of the A 2 singularity. The width of the gap of the dis¬ 
continuous jump of P decreases as c is increased. As c 
increases further, the situation qualitatively changes. The 
minimum of p{B, c) at Bd becomes unstable and eventu¬ 
ally disappears. It is clear from eq. ([T2j) that this happens 
at Cg = 1/5. At this point, the A 2 singularity disappears 
and the transition becomes continuous. This is the signal 
of the As singularity. The MCT predicts that the criti¬ 
cal behavior of the A 3 singularity is distinct from that of 
the A 2 singularity [3j|6j[7l|9] . Qualitatively the same trend 
is observed for different sets of k and / (see the inset of 
Fig. [2]). 

p{c,p) evaluated by numerically solving eqs. (| 8 |) — ([TT]) is 
plotted in Fig. [3] for several c’s. The behavior of p{c,p) is 
qualitatively the same as that of the nonergodic parameter 
evaluated from the MCT across the A 3 singularity point. 
Note that similar A 3 singularities have been analyzed for 
binary and ternary mixtures of the FAM on the Bethe 
lattices and random graphs HSU]. 
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Fig. 4: The (dynamic) phase diagram of the randomly pinned 
FAM: The thick solid line is the nonergodic transition line 
Pd{c). The dashed line is a “spinodal” line at which the larger 
branches of B disappear. The filled circle at the end of the 
thick solid line at (c = 1/5, p = 5/6) is the terminal critical 
point of the higher order singularity. The arrows represent the 
directions of vectors e± and ey (see the text). 


The nonergodic transition line as a function of c, Pd{c), 
is shown in Fig. 0] in the thick solid line. Recall that there 
exist three solutions of B of eq. m if p is slightly above 
Pd and the smallest B is the real solution. As it is ob¬ 
served in Fig. [21 if c is slightly smaller than Cg, there is 
always a finite p {> Pd) above which the larger solutions 
of B disappear. In Fig. 01 this “spinodal” line is shown 
in dashed line. This “phase-diagram” demonstrates that 
the transition point pd decreases (or Td increases) as c in¬ 
creases and eventually terminates at a critical point. The 
critical point can be evaluated analytically as = 1/5 
and Ps = Pt^(cs) = 5/6. At c > c^, there is no distinc¬ 
tion between the ergodic (fluid) and nonergodic (glassy) 
phases and (j) continuously increases as p increases. The 
behavior of Pd{c) in Fig.0]is qualitatively the same as that 
evaluated for the randomly pinned p-spin mean-field spin 
glass model |6ir7| . We present the phase diagram for /c = 3 
and / = 2 but the results remain unchanged for other sets 
of k and /, as long as /c > / > 1 . 


Higher Order Singularity. — If c < c^, the transi¬ 
tion across the thick solid line in Fig. 0] is discontinuous 
and the order parameters, (j) (or 5), behave as (/) — (pd (or 
B — Bd) oc \pd — [IS]. At c = Cg = 1/5, however, the 

properties of the transition qualitatively change. In the 
following, we show that the critical behavior around this 
terminal critical point is indeed characterized as the A 3 
singularity [3l[9l [2^[23] . 

We first introduce a function defined by 


Q{c,p,B) = l-p+{l-c)p{B^ + B)}-B. (13) 


From eq. m, Q = 0 always holds. We shall expand this 


function around the transition point as 


Q{c,p,B)=^5B + 

+ 

op 


^ A 


1 d^Q 
d^Q 

dpdB^ 


SB'^ 


dc 


Sc 


5p] SB 


(14) 


where SB = B — Bd and Sc = c — Cd. The transition 
point is called the A 2 singularity point if dQ/dB = 0 and 
d^QfdB‘^ 7 ^ 0, and the A 3 singularity point if dQjdB = 
d^QjdBA = 0 and d^QjdBA ^ 0 [91 I221I23] . It is obvious 
that on the transition line, Pd{c), of Fig. 0] for small c, 
dQjdB = 0 and d^Q/dBA ^ 0 are satisfied and that they 
are of the A 2 type. However, exactly on the terminal criti¬ 
cal point, (cs,Ps,Hs) = (1/5, 5/6,1/2), vanishes 

and this point is categorized as the A 3 singularity point. 
In the vicinity of the terminal critical point, eq. (pA]) can 
be written as 


where we defined 5c = {cs — cjjcs and = (ps — p)lps- 
Note that we have two control parameters, s = {Sc^Sp). 
In order to investigate the critical behavior, it is conve¬ 
nient to transform the set of control parameters parallel 
and perpendicular to the transition line. Since the first 
two terms of eq. m are written in the form of the prod¬ 
uct of e and a constant vector (1/12,1/2), it is natural 
to introduce the unit vectors directing to and perpen¬ 
dicular to this constant vector by ej_ = -^(1,6) and 
ey = (—6,1) (see arrows in Fig. 0|), so that s can be 

written as 5 = e±e± -h5||e||. With this expression, eq. m 
can be rewritten as 

iSB^ , 3<5R(25£x-2£||) 

3 12 16v^ ’ ^ ^ 

and one immediately finds SB ^ if 7 ^ 0. On the 

other hand, if 5 _l = 0, or if the critical point is approached 
from the direction parallel to ey, one finds SB ^ (—5y)^/^. 


Correlation lengths. — It is established, in the con¬ 
text of the mean-field scenario of the glass transition, 
that the cooperative dynamics are associated with diverg¬ 
ing correlation lengths near the dynamic glass transition 
point. These lengths are predicted by the inhomogeneous 
version of the MCT (IMCT) both for the A 2 and A 3 sin¬ 
gularity points [9l [24ll26] . In this section, we discuss the 
correlation lengths of the randomly pinned FAM, using a 
point-to-set function, which is a quantity introduced to 
analyze the spatial correlations in glassy systems [lilEH]. 

The local structure of the Bethe lattice is Cayley-tree¬ 
like as shown in Fig. [H We freeze the degrees of freedom of 
the spins at the outermost branch do and set Bq = 0. We 
consider how much the effect of the frozen boundary pen¬ 
etrates down to the inner branches. The value of Bi down 
from the 0 -th node is written by the recursive equation as 

Ri+i = l-p+{l-c)p {Rf + 3R|(1 -Bi)]. (17) 
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Fig. 5: Dependence of Bi on the depth of the inner branch 
for c = 0 and for various p’s. We define {pd —p)lpd = 2“^^ and 
the lines with circles correspond to n = 1, 2, • • • , 10 from left 
to right. The dashed horizontal line represents the threshold 
value 1/4. 

In Fig.O the numerical solution for Bi of eQ. (fT7|) for c = 0 
is shown as a function of i for several values of p’s (< pd) 
in ergodic states. Bi for i > 0 remains very small near 
the boundary but it relaxes to 1 far away from the bound¬ 
ary (i 1 ). As p approaches to pd from below, the dis¬ 
tance over which the effect of the boundary penetrates 
grows. We shall define the correlation length ^ as the i- 
th point at which Bi exceeds the threshold whose value 
is set to Bd, the solution of eg. (ITT]) at p = p^^ [29]. Nu¬ 
merical results for ^ for c = 0 , 0 . 1 , and 0.2 are shown 
as a function of the distance from the transition point 
^ = {Pd —p)IPd in Fig- El It is clearly seen that ^’s behave 
as \pd — for c = 0 and 0.1 [29]. But for c = 1/5, 

one observes ^ oc \pd — This asymptotic behavior 

of ^ can be analytically explained by the linear analysis 
of eq. (El) around pd- Let us define the distance of Bi 
from the plateau value Bd by SBi = Bi — Bd- We assume 
that the variation of Bi is small and can be represented 
using a large number I as SBi_^i = (1 — l~^)SBi. Since 
6Bi ^ (1 — ^ I can be identified as p Sub¬ 

stituting this expression back to the recursive equation, 
eq. El , we arrive at 

C^SBi = Qic,p,Bi), (18) 

where Q{c,p^Bi) is defined by eq.([T3|). Expanding 
Q{c,p^Bi) in terms of SBi around Bd like eq. flH]) . it is 
straightforward to show that ^ behaves as ^ ex 5 “^/^ 
when c < 1/5. In the vicinity of the terminal critical 
point (cs,Ps) = (1/5, 5/6), the behavior of ^ qualitatively 
changes. If 5 _l p 0, we have ^ (x thus reproducing 

the results of Fig.O If 5 _l = 0, one finds ^ (x 5 ^^. The di¬ 
rection dependence of the growth of ^ is shown in Fig.Qto- 
gether with the direct numerical results. It is known that 
the exponents for the length scales of the model systems 
on the Bethe lattices are twice as large as the correspond¬ 
ing values of the mean-field theory [29H3T] . Therefore, 
the exponent of the diverging length ^ ^ s~^ obtained 
above should correspond to the value of the mean-field 
theory of u = 1/4 instead of 1/2 for c < 1/5. Likewise, 


10 ^ 

10 ^ 

10 ^ 

10 ° 

10-6 10-4 ^ 10-2 10 ° 

Fig. 6: ^ as a function oi e — {pd — p)/Pd for several c’s. Circles 
are for c = 0, triangles for c = 0.1, and squares are for c = 0.2. 
Thin dotted and dashed lines represent and 

respectively. 

near the terminal critical point, the exponents should be 
= 1/3 for 5 _l p 0 and z/y = 1/2 for 5 _l = 0, respectively. 
These results are equivalent with those obtained from the 

IMCT [HES]. 

Conclusion. — In this letter, the stationary proper¬ 
ties of the randomly pinned FAM on the Bethe lattice have 
been discussed. The order parameter, 0, or the probabil¬ 
ity that the spin flips, 5, exhibits analogous behavior as 
the nonergodic parameter of the mean-field models of the 
glass transition in the presence of the randomly pinned 
spins/atoms at their dynamic transition lines [5]-[7] . As the 
fraction of the pinned spins, c, increases, the A 2 singular 
point Pd{c), at which the order parameter discontinuously 
jumps, decreases (or Td{c) increases) and terminates at a 
larger but finite c = Cg, where the transition point be¬ 
comes the higher order singular point of type A 3 , which 
is fully consistent with the dynamic transition lines pre¬ 
dicted by the MCT for randomly pinned systems. The 
concept of the point-to-set correlation has been employed 
to evaluate the correlation length, around this transition 
line, Pd{c) [27l[3T]. Aside from the trivial factor of 2 , the 
exponent, z/, of the correlation lengths, ^ ^ 5 “^, is consis¬ 
tent with that obtained from the inhomogeneous version 
of the MCT [9] . The exponent at the A 3 singularity point 
is = 1/3 and z/y = 1/2, depending on the directions 
to approach the critical point. This behavior is akin to 
that of the ferromagnetic transition in the Landau theory 
(z/ = 1/2 when the magnetic field h = 0 and T ^ Tc and 
z/ = 1/3 when T = Tc and h ^ 0 |9]). Note that the 
parallel direction is more special than the perpendicular 
direction, because the critical behavior along the former 
is observed only when = 0, whereas the later is generic 
in a sense that the critical behavior is observed even if 
£± ^ 0. We should emphasize, however, that the analysis 
of the order parameter as studied here can not establish 
the universality class of the transition. It is known that 
the transition of the mean-field random-field Ising model 
(RFIM) is also characterized by the same exponents as 
those of the ferromagnetic transition in the Landau the¬ 
ory but its universality class is distinct [32]. The analysis 


c = 0.2 

‘i-i- 

c = 0.0 ^ 
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Fig. 7: ^ near the terminal critical point (cs = 1/5, ps = 
5/6) approached from the parallel (J_) and perpendicular (||) 
directions. Numerical results are shown in filled symbols. Thin 
dashed lines represent and |£||| , respectively. 


of the fluctuations such as the higher order moments is 
necessary to characterize the nature of the transition. In¬ 
deed, it has been demonstrated that from the dynamic 
scaling of the model discussed here at c = 0 , the model 
belongs to the RFIM universality class [33]. Besides, the 
randomly pinned glass model is also argued to belong to 
the same universality class [ 8 ]. 

In this letter, we only analyzed the stationary proper¬ 
ties of the randomly pinned FAM. It is known that the 
slow relaxation dynamics of the FAM on the Bethe lat¬ 
tice for c = 0 is qualitatively similar to that of the MCT 
for the bulk supercooled liquid or mean-field spin glass 
models [16]. On the other hand, the MCT and IMCT 
predict that the slow dynamics near the A 3 singularity 
point is distinct from that at the A 2 point, and the sys¬ 
tem exhibits anomalous dynamics such as the logarithmic 
relaxation and milder growth of the correlated fluctuations 
[319]. Studies in this direction are left for future works. 
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